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ABSTRACT 

The second multisite photometric campaign devoted to v Eri is reported. The cam- 
paign, carried out from 11 Sept. 2003 to 16 Feb. 2004, was very nearly a replica of the 
first, 2002-3 one: the five telescopes and photometers we used were the same as those 
in the first campaign, the comparison stars and observing procedure were identical, 
and the numbers and time base-lines of the data were comparable. 

For v Eri, analysis of the new data adds four independent frequencies to the nine 
derived previously from the 2002-3 data, three in the range from 7.20 to 7.93 d _1 , 
and a low one, equal to 0.614 d _1 . Combining the new and the old data results in 
two further independent frequencies, equal to 6.7322 and 6.2236 d" 1 . Altogether, the 
oscillation spectrum is shown to consist of 12 high frequencies and two low ones. The 
latter have u amplitudes about twice as large as the v and y amplitudes, a signature 
of high radial-order g modes. Thus, the suggestion that v Eri is both a (3 Cephei and 
an SPB star, put forward on the basis of the first campaign's data, is confirmed. 

Nine of the 12 high frequencies form three triplets, of which two are new. The 
triplets represent rotationally split 1=1 modes, although in case of the smallest- 
amplitude one this may be questioned. Mean separations and asymmetries of the 
triplets are derived with accuracy sufficient for meaningful comparison with models. 

The first comparison star, fj, Eri, is shown to be an SPB variable with an oscillation 
spectrum consisting of six frequencies, three of which are equidistant in period. The 
star is also found to be an eclipsing variable. The eclipse is a transit, probably total, 
the secondary is fainter than the primary by several magnitudes, and the system is 
widely detached. 

The second comparison star, £ Eri, is confirmed to be a S Scuti variable. To the 
frequency of 10.8742 d -1 seen already in the first campaign's data, another one, equal 
to 17.2524 d" 1 , is added. 

Key words: techniques: photometric - stars: early-type - stars: individual: v Eridani 
- stars: individual: fi Eridani - stars: individual: £ Eridani - stars: oscillations - stars: 
eclipsing 
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1 INTRODUCTION 

The first multisite photometric and spectrographic cam- 
paign devoted to the (3 Cephei star v Eridani was carried 
out between October 2002 and February 2003. More than 
600 h of differential uvyV photometry on 148 nights and 
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more than 2000 high-resolution spectra were obtained. A 
fre quency analysis of t he photometric data was reported 
bv lHandler et all (|2004L hereafter Paper I), while the spec- 
trographs t i me-se ries and their analysis were presented by 
lAerts et alj i2004l hereafter Paper II). An extended fre- 
quency analysis and m ode identification was provided by 
iDe Ridder et alJ <2004 hereafter Paper III). 

Seismic modelling of the oscillation spec- 
trum of v Eri has been u nderta ken by 
Pajny_atn^ddi J Handler fc Dziembowskil i2004) and 
Ausseloos et alJ <2004T ~ 

In Paper I, the light-variation of v Eri was shown to 
consist of 23 sinusoidal terms. These included 8 independent 
ones with frequencies, fi, i — 1..8, spanning the range from 
5.6 to 7.9 d _1 , 14 high-frequency combination terms, and a 
term with the low frequency /a = 0.432 d _1 . 

The f our h ighest-amplitude terms, discovered long ago 
bv ISaitol dl97d) . consist of a singlet with frequency f± = 
5.763 d _1 , and a triplet very nearly equidistant in frequency 
(/ 3 = 5.620, fi = 5.637 and f 2 = 5.654 d _1 ). Before the 
campaign it was known that the singlet wa s a radial mode 
llCugier. Dziembowski fc Pamya tnvkhll l994l). and surmised 
that it was the fundamental ([Dziembowski fc .Terzvkiewiczl 
2003), but the true nat ure of the triplet was unclear. From 
multic olour photometry. iHevnderickx. Waelkens fc Smeversl 
(1994) derived I = 1 for f<t, and suggested that the triplet 
wa s a rotationally split dipole mode. This was challenged 
bv lAerts. Waelkens fc de Pauwl (1 19941) who — from an anal- 
ysis of line-profile variations — identified the f% term with an 
axisymmetric mode. However, both results are questionable 
because in neither case the triplet had been resolved. The 
matter was settled in Paper III: the wavelength dependence 
of the uvy amplitudes of the triplet t erms implies I = 1 for 
all of them. IPamvatnvkh et aD ll2004l) showed then that the 
triplet was a g\ mode. 

A frequency triplet /_ < fo < f+ can be characterized 
by its mean separation 

S = 0.5(/+ - /_), (1) 
and asymmetry 

A = f-+f+-2f . (2) 

In case of a rotationally split triplet, S is determined by the 
angular rotation rate of the star, fi, while A is sensitive to 
effects of higher order in fi, as well as to effects of a magnetic 
field. If these effects are negligible, A = 0. 

For the i — 3, 4, 2 triplet, both parameters were 
derived before the campa i gn fr om archival data by 
iDziembowski fc Jerzvkiewiczl ll2003l) . The value of asym- 
metry they obtained, A = (-7.1 ± 0.3) x 10" 4 d" 1 , was 
unexpectedly large for the s mall rotation rate implied by 
the tr iplet's S = 0.0168 d -1 . IDziembowski fc Jerzvkiewiczl 
(2003) showed that this problem may be solved by postulat- 
ing the existence of a 5 to 10 kG magnetic field in the outer 
envelope of the star. From the campaign data (see Paper 
III), one gets A = (— 4.9 ± 1.1) x 10~ 4 d~ x , a value which is 
not in serious conflict with the observed S. Unfortunately, 
the large standard deviation of this result makes it useless. 
A longer time base-line than that of the 2002-3 campaign 
would be needed to obtain a more reliable value of A and 
thus decide whether invoking magnetic field were necessary. 



This was one motivation for undertaking the sequel cam- 
paign. 

In Paper III, the spherical harmonic degree of the i = 5, 
6 and 7 terms was found to be I = 1, but an attempt to de- 
rive I for the low-frequency i = A term (referred to as the 
vio term in that paper) was unsuccessful because of its small 
uvy amplitudes and the poor resolution of diagno stic dia- 
grams at low frequencies. IPamvatnvkh et al.l ll2004f) showed 
that the i = 5 term is a p2 mode, while the i = 6 one is a 
pi mode. Moreover, they suggested that the low-frequency 
term is an i = 1, m = — 1, <?i6 mode. 

According to IPamvatnvkh et al.l (I2004T) . only the i = 1, 
2, 3 and 4 modes (i.e., the radial fundamental mode and the 
I = 1, gi triplet) are unstable in standard models. Modes 
with / > 6 d _ (i.e., the i = 5, 6 and 7 ones) and the 
low-frequency mode are stable. IPamvatnvkh et alJ (120041) 
demonstrate, however, that a fourfold overabundance of Fe 
in the driving zone would account for excitation of the high- 
frequency modes, and would make the low- frequency mode 
marginally unstable. 

It has been noted in Paper I that if the low-frequency 
term were indeed a high-order g mode, v Eri would be both 
a Cephei variable and a slowly pulsating (SPB) star. Un- 
fortunately, /a differs from the sixth order combination fre- 
quency 3/i —3/3 by less than 0.003 d" 1 . This number is much 
larger than the formal error of /a, but is smaller than half 
the frequency resolution of the campaign data. Thus, the 
possibility that /a is the combination frequency — although 
rather unlikely — cannot be rejected. Again, a longer time 
base-line would help to settle the issue. 

In addition to extending the time base-line, the sequel 
campaign was expected to double the number of data points 
and thus lower the detection threshold so that modes hav- 
ing amplitudes too low to be seen in the 2002-3 frequency 
spectra could be discovered from the combined data. 

Both comparison stars used in the 2002-3 photometric 
observations turned out to be variable. For the first, fi Eri 
(HD 30211, B5 IV, V = 4.00), the analysis carried out in Pa- 
per I revealed a dominant frequency f{ = 0.6164 d _1 . (From 
now on we shall use a prime to denote frequencies of /j, Eri.) 
Prewhitening with this frequency resulted in an amplitude 
spectrum with a very strong 1/f component, indicating a 
complex variation. Since the star is a spe ctroscopic binary 
with an orbital period P orb = 7.35890 d llHilll 1 1968ft . tests 
have been made to detect a non-sinusoidal signal with this 
period. Unfortunately, they were inconclusive. 

Taking into account the star's position in the HR. dia- 
gram, the frequency f[ of the dominant variation, and the 
fact that the u amplitude was about a factor of two greater 
than the v and y amplitudes, it was concluded in Paper I 
that /1 Eri is probably an SPB star. However, it was also 
noted that instead of pulsation, a rotational modulation 
could be the cause of the dominant variation. 

In Paper I, the frequency analysis of de-trended differ- 
ential magnitudes of \x Eri revealed a small-amplitude vari- 
ation with a frequency / x = 10.873 d _1 . It was suggested 
that the second comparison star, £ Eri (HD 27861, A2V, 
V = 5.17), may be responsible. 

The present paper reports the sequel photometric cam- 
paign, carried out from 11 Sept. 2003 to 16 Feb. 2004, and 
an analysis of the data of both campaigns. The 2003-4 ob- 
servations and reductions are described in Sect. 2. Sect. 3 
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contains an account of the frequency analysis of the new 
data and a comparison of the results with those of Paper I. 
Low frequencies in the variation of \jl and v Eri from the uvy 
data of the first campaign are re-examined in Sect. 4. Sect. 
5 is devoted to frequency analysis of the combined, 2002-4 
data. Finally, Sect. 6 provides a summary with an emphasis 
on clues for asteroseismology of the three stars, v, /i and £ 
Eri. 



2 OBSERVATIONS AND REDUCTIONS 

The observations were carried out with five telescopes on 
four continents. An observing log is presented in Table 1. 
Comparison of this table with Table 1 of Paper I shows 
that the new data are almost as extensive as the old ones. 
However, the 2003-4 y data are less numerous than the v 
and u ones (see below). The time base-lines of the 2002-3 
and 2003-4 sets are very nearly the same; they amount to 
157.9 and 158.5 d for the old and new data, respectively. 

The five telescopes and photometers were the same as 
those in the 2002-3 campaign. Thus, single-channel photo- 
electric photometers were used at all sites but Sierra Nevada 
Observatory (OSN), where a simultaneous uvby photome- 
ter was used. At OSN, the observations were obtained with 
all four Stromgren filters, at Fairborn, Lowell and Siding 
Spring, with u, v and y. At SAAO, the y filter used in the 
2002-3 campaign has deteriorated to such a degree that it 
had to be discarded. Since no replacement was available, the 
data were taken with two filters, Stromgren u and v, except 
that on his first two nights AP used Johnson filters B and 
V. 

The comparison stars and observing procedures were 
the same as in the first campaign. 

The reductions, carried out separately for each of the 
three wavelength bands, u, v, and y or V, consisted in 
(1) computing heliocentric JD numbers for the mean epochs 
of observations, (2) computing the air mass for each obser- 
vation, (3) correcting instrumental magnitudes of v, fi and 
£ Eri for atmospheric extinction with first-order extinction 
coefficients derived from the instrumental magnitudes of £ 
Eri by means of Bouguer plots, (4) forming differential mag- 
nitudes l u Eri — £ Eri' and '/x Eri — £ Eri,' and (5) setting 
the mean light-levels of the differential magnitudes from dif- 
ferent telescope-filter combinations to the same values. In 
step 3, second-order extinction corrections were not applied 
because no colour-dependent extinction effects could be de- 
tected in the uncorrected differential magnitudes (but see 
the last paragraph of Sect. 3.1). In step 4, the magnitudes 
of £ Eri were interpolated on the epoch of observation of v or 
/x Eri. In step 5, the mean-light levels for each telescope-filter 
combination were derived using residuals from least-squares 
solutions with the four highest-amplitude terms (see the In- 
troduction). 



3 FREQUENCY ANALYSIS OF THE NEW 
DATA 

The analysis was carried out in essentially the same way as 
in Paper I. That is, frequencies were identified from peri- 
odograms, one at a time. Before each run, the data were 



prewhitened with all frequencies found in previous runs. Af- 
ter several runs, the frequencies were refined by means of a 
nonlinear least-squares fit using the values of independent 
frequencies read off the periodograms as starting values. 
The frequencies of the combination terms were computed 
from the independent frequencies. Thus, the unknowns in 
the normal equations were the corrections to the indepen- 
dent frequencies, to the mean magnitude, (Am), and to the 
amplitudes, Ai, and phases, <j>i, appearing in the following 
expression: 

JV 

Am = (Am) + ^ A 1 sm[2n f t (t - to) + (3) 

i=l 

where Am is the differential magnitude in u, v, or yV, N 
is the number of all frequencies, fi, the combination terms 
included, and to is an arbitrary initial epoch. 

The difference with respect to the 2002-3 analysis con- 
sisted in using different programs: the PERIOD 98 package 
(see Paper I) was replaced with programs that have been in 
use by MJ since 1975 (see, e.g.. |Jerzvkiewic^ll978l) . Thus, 
by "periodogram" we now mean a power spectrum, and not 
an amplitude spectrum as in Paper I; by "power" we mean 
1 — a 2 (f)/<j 2 , where & 2 (f) is the variance of a least-squares 
fit of a sine-curve of frequency / to the data, and a 2 is the 
variance of the same data. However, for the purpose of esti- 
mating signal-to-noise ratios (see below) we also computed 
amplitude spectra. 

3.1 The programme star 

The data used for analysis were the differential magnitudes 
l v Eri — £ Eri.' Power spectra were computed independently 
from the u and v data. No power spectra were computed 
from the less numerous yV data, but nonlinear least-squares 
solutions were carried out for all three bands, with the start- 
ing values of the independent frequencies for y (henceforth 
we shall use "y" instead of "yV") taken from v. The OSN b 
and SAAO B data were not used. 

In the u and v power spectra, 14 independent and 15 
combination frequencies could be clearly seen. They are 
identified in the first column of Table 2. The numbers in 
the table are from nonlinear least-squares solutions. The 
values of the independent frequencies and their standard 
deviations, computed as straight means from the separate 
solutions for u, v and y, are given in column 2 above the 
horizontal line. The combination frequencies, listed below 
the line, were computed from the independent frequencies 
according to ID in the first column; their standard devia- 
tions were computed from the standard deviations of the in- 
dependent frequencies assuming rms propagation of errors. 
The amplitudes, A u , A v and A y , with the standard devi- 
ations, given in columns 3, 4 and 5, respectively, are from 
the independent solutions for u, v and y. The last column 
lists the v-amplitude signal-to-noise ratio, S/N, defined as 
in Paper I, except that in Paper I the mean noise level was 
estimated in 5 d _1 intervals, while now we adopted 0.1 d _1 
intervals for frequencies lower than 3 d _1 , and 1 d _1 for 
higher frequencies. In all cases, the amplitude spectra of the 
data prewhitened with the 29 frequencies of Table 2 were 
used for estimating the mean noise level. 

In Paper I, a peak in the amplitude spectrum was con- 
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Table 1. Log of the photometric measurements of v Eri in 2003-4. Observatories are listed in the order of their 
geographic longitude. 



Observatory 


Longitude 


Latitude 


Telescope 


Amount of data 


Observer (s) 










Nights 


h 




Sierra Nevada Observatory 


-3°23' 


+37°04' 


0.9-m 


4 


10.64 


ER 


Fairborn Observatory 


-110° 42' 


+31°23' 


0.75-m APT 


51 


211.09 




Lowell Observatory 


-111°40' 


+35°12' 


0.5-m 


25 


87.10 


MJ 


Siding Spring Observatory 


+149°04' 


-31°16' 


0.6-m 


26 


79.77 


RRS 


South African Astronomical Observatory 


+20° 49' 


-32°22' 


0.5-m 


17 


92.13 


AP 


South African Astronomical Observatory 


+20° 49' 


-32°22' 


0.5-m 


19 


48.52 


RM, TM, PT 


Total 








142 


529.25 





Table 2. Frequencies and amplitudes in the differential magnitudes V Eri — £ Eri' from the 
2003-4 data. Independent frequencies are listed in the upper part of the table. The combination 
frequencies are listed below the horizontal line. In both cases the frequencies are ordered according 
to decreasing v amplitude, A v . The last column contains the ^-amplitude signal-to-noise ratio. 
Frequencies / x and f y are due to f Eri. 



ID Frequency [d 1 ] A u [mmag] A v [mmag] A y [mmag] S/N 





h 


5.763256 ± 





.000012 


72. 


3 


± 


0. 


19 


40. 


.6 


± 


0. 


14 


36 


.7 


± 


0. 


15 


165.0 




k 


5.653897 ± 





.000020 


38. 


6 


± 


0. 


19 


27 


.2 


± 


0. 


14 


25. 


A 


± 


0. 


15 


110.6 




h 


5.619979 ± 





.000021 


35. 


5 


± 


0. 


19 


25 


.0 


± 


0. 


14 


24. 


.0 


± 


0. 


15 


101.6 




h 


5.637215 ± 





.000025 


31. 





± 


0. 


18 


21. 


.8 


± 


0. 


13 


20 


A 


± 


0. 


16 


88.6 




h 


7.89859 


± 





.00022 


3 


.7 


± 


0. 


19 


2 


.7 


± 


0. 


14 


2. 


,5 


± 


0. 


16 


11.0 




h 


6.26225 


± 





.00025 


3 


.1 


± 


0. 


19 


2 


.2 


± 


0. 


14 


2. 


,2 


± 


0. 


16 


8.8 




/a 


0.43257 


± 


0. 


.00028 


3. 


.2 


± 


0. 


19 


2 


.0 


± 


0. 


14 


1. 


.7 


± 


0. 


16 


5.0 




h 


6.24468 


± 





.00034 


2. 


.3 


± 


0. 


18 


1 


.6 


± 


0. 


13 


1. 


.6 


± 


0. 


15 


6.4 




/b 


0.61411 


± 





.00033 


3 


.4 


± 


0. 


18 


1 


.5 


± 


0. 


13 


1. 


.6 


± 


0. 


15 


4.3 




/x 


10.8743 


± 


0. 


.0006 


1. 


.2 


± 


0. 


18 


1 


.4 


± 


0. 


13 


0. 


.9 


± 


0. 


15 


5.7 




/io 


7.9296 


± 


0. 


.0005 


1. 


.6 


± 


0. 


19 


1 


.2 


± 


0. 


11 


1 


.2 


± 


0. 


16 


5.3 




h 


7.9132 


± 





.0005 


1. 


.8 


± 


0. 


18 


1 


.1 


± 


0. 


13 


1. 


.3 


± 


0. 


16 


4.9 




h 


7.2006 


± 





.0005 


1. 


.4 


± 


0. 


18 


1 


.0 


± 


0. 


13 


1. 


.1 


± 


0. 


15 


4.4 




fy 


17.2534 


± 


0. 


.0006 


1. 





± 


0. 


19 





.9 


± 


0. 


13 


0. 


.9 


± 


0. 


15 


1.6 


h 


+ /2 


11.417153 


± 





.000023 


12. 


3 


± 


0. 


19 


8 


.8 


± 


0. 


14 


8. 


.4 


± 


0. 


15 


34.0 


h 


+ h 


11.383235 


± 





.000024 


10. 


.8 


± 


0. 


19 


7 


.6 


± 


0. 


11 


7 


.1 


± 


0. 


15 


29.3 


h 


+ /4 


11.400471 


± 





.000028 


9. 


.9 


± 


0. 


19 


7 


.0 


± 


0. 


11 


6. 


.6 


± 


0. 


16 


27.0 




2/i 


11.526512 


± 





.000017 


4. 


.5 


± 


0. 


18 


3 


.1 


± 


0. 


13 


3 


.1 


± 


0. 


15 


12.0 


+ /2 


+ /3 


17.037132 


± 





.000031 


3 


.3 


± 


0. 


18 


2 


.2 


± 


0. 


13 


2. 


.0 


± 


0. 


15 


11.3 


h 


-h 


0.109359 


± 





.000023 


2 


.7 


± 


0. 


19 


1 


.6 


± 


0. 


14 


1. 


.6 


± 


0. 


16 


4.4 


2/i 


+ /4 


17.163727 


± 


0. 


.000030 


1. 


.9 


± 


0. 


IS 


1 


.4 


± 


0. 


13 


1. 


.3 


± 


0. 


15 


7.2 


2/i 


+ /2 


17.180409 


± 





.000026 


1. 


.9 


± 


0. 


18 


1 


.4 


± 


0. 


13 


1. 


,3 


± 


0. 


15 


7.2 


h 


+ h 


11.273876 


± 





.000029 


2. 


.4 


± 


0. 


19 


1 


.3 


± 


0. 


14 


1. 


,1 


± 


0. 


15 


5.0 


2/i 


+ h 


17.146491 


± 


0. 


.000027 


1. 


.4 


± 


0. 


18 


1 


.2 


± 


0. 


13 


0. 


.8 


± 


0. 


15 


6.2 




2/2 


11.307794 


± 


0. 


.000028 


1. 


.1 


± 


0. 


19 


1 


.1 


± 


0. 


14 


0. 


.6 


± 


0. 


16 


4.2 


h 


-h 


0.126041 


± 





.000028 


1. 


.9 


± 


0. 


18 


1 


.1 


± 


0. 


13 


1. 


.2 


± 


0. 


15 


3.0 


+ h 


+ /4 


17.020450 


± 





.000035 


1. 


.3 


± 


0. 


18 


1 


.1 


± 


0. 


13 


0. 


.9 


± 


0. 


15 


5.6 


+ h 


+ /3 


22.800388 


± 


0. 


.000034 


1. 


.3 


± 


0. 


18 


1 


.1 


± 


0. 


13 


0. 


.8 


± 


0. 


15 


6.1 


+ h 


-h 


5.797174 


± 





.000031 


1. 


.5 


± 


0. 


19 





.9 


± 


0. 


11 


1 


.3 


± 


0. 


15 


3.7 



sidered to be significant if its signal-to-noise ratio exceeded 
4 for an independent term or 3.5 for a combination term. 
As can be seen from Table 2, this condition is met by all 
terms identified from the 2003-4 data except the differential 
combination term /i — f$ for which S/N = 3.0. This term 
may be spurious. 

The standard deviations in Table 2 (and in the tables 
that follow), referred to as "formal" in the remainder this 
subsection and in the first three paragraphs of the next sub- 



section, will be underestimated if the errors of differential 
magnitudes are correlated in time. For the case at hand, 
i.e., of fitting a sinusoi d to time series data, the pro blem 
has been discussed by ISchwar zenberg^ Czernvl <199lfl and 
iMonteomerv fc O'DonoghueTi^OOr) . These authors consider 
the factor, D, by which the formal standard deviations of a 
frequency, amplitude and phase should be multiplied in or- 
der to get correct values. (D is the same for frequencies, 
amplitudes and phases because the ratio of the formal stan- 



dard deviations of any two of these quantities is determined 
by error-free numbers such as the epochs of observations and 
the starting values of the frequencies.) The factor depends on 
an estimate of the number of consecutive data points which 
are correlated. Unfortunately, this estimate is not easy to 
obtain, especially for time series such as ours, consisting of 
data from many nights and several observatories. We shall 
return to this point in the next subsection. 

Figure shows the power spectra of the u and v differ- 
ential magnitudes l v Eri — £ Eri' prewhitened with the 29 
frequencies of Table 2. In both panels of the figure, the high- 
est peaks are seen at frequencies lower than about 5 d~ .In 
the u spectrum (lower panel), the highest peak occurs at 
0.314 d , while in the v spectrum, the highest peak is the 
one at 1.498 d _1 . The u amplitude at 0.314 d" 1 is equal to 
2.2 mmag, while the v amplitude at 1.498 d _ , to 1.1 mmag. 
The u- and D-amplitude S/N amount to 2.7 in both cases. 
We conclude that these frequencies are probably spurious. 

For frequencies higher than 5 d _1 , the highest peak in 
the u spectrum (lower panel) occurs at 5.015 d _1 . For this 
peak the u-amplitude S/N is equal to 5.0, so that the signifi- 
cance condition mentioned earlier in this section is satisfied. 
However, there is no power at 5.015 d" 1 in the v spectrum 
(upper panel); the peak closest to 5.015 d _1 , at 5.128 d _1 , 
is probably noise because it has S/N equal to 3.2. This, and 
the fact that 5.015 d _1 is very nearly equal to 5 cycles per 
sidereal day (sd _1 ), suggest that the peak is due to colour 
extinction in the u band which we neglected (see Sect. 2). 
Indeed, the differential colour-extinction correction contains 
a term equal to k"XAC, where k" is the scond-order extinc- 
tion coefficient, X is the air-mass, and AC is the differen- 
tial colour-index. Because of the second factor, X, the term 
causes a parabola-shaped variation symmetrical around the 
time of meridian passage. For a single observatory, this will 
produce a signal with frequency equal to n sd~\ where n 
is a small whole number equal to the number of sinewaves 
that best fit the variations on successive nights. In our case 
this number is apparently equal to 5. For multi-site data, 
phase smearing may occur, tending to wash out the signal. 
However, our time-series is dominated by the data from Fair- 
born and Lowell, two observatories lying on nearly the same 
meridian and therefore introducing negligible phase shift. 
In addition, the greatest common divisor of the longitude 
differences between Fairborn and Lowell on one hand, and 
SAAO and Siding Spring on the other is close to 60°, making 
the 5 sd _1 signals from these observatories to approximately 
agree in phase. We conclude that the 5.015 d _1 peak in the 
u spectrum in Fig. Q is not due to an intrinsic variation of 
v Eri. 

3.2 Comparison of the 2003-4 and 2002-3 results 
for the programme star 

Of the 14 independent terms in Table 2, nine appear in Table 

2 of Paper I. The differences between the 2003-4 and 2002- 

3 values of their frequencies and amplitudes are listed in 
Table 3. The standard deviations given in the table were 
computed from the standard deviations of the 2003-4 and 
2002-3 least-squares solutions assuming rms propagation of 
errors. Thus, they are formal in the sense defined in the 
previous subsection. 

We shall now return to the factor D by which the for- 
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Figure 1. Power spectra of the u (lower panel) and v (upper) 
differential magnitudes 'v Eri — £ Eri' prewhitened with the 29 
frequencies of Table 2. 



mal standard deviations are underestimated. Let us con- 
sider the nine independent frequencies common to 2002-3 
and 2003-4. As can be seen from Table 3, the moduli of 
the frequency differences, |A/|, are of the same order of 
magnitude as the formal standard deviations of A/, <TAf. In 
four cases |A/| < <TAf, in three cases (TAf < |A/| < 1o Af , 
and in two cases the differences exceed 2<7Af • Taken at face 
value, these inequalities would indicate that four or five of 
the nine frequencies changed from 2002-3 to 2003-4 by lo"Af 
or more, while /3, /s and fs changed by 2<7Af or more. How- 
ever, frequency changes in (3 Cephei stars have time scales 
much longer than 1 year. For the fou r first frequenc i es of v 
Eri, this is shown to be the case by lHandler et af] (2005). 
Although the existence of long-term variations does not ex- 
clude the possibility of year-to-year ones, let us assume that 
these four frequencies were strictly constant from 2002 to 
2004. If the assumption were false, the value of D we derive 
in the next paragraph will be too large. 

For / = const, where / is any of the four frequencies, 
the modulus of Af can be thought of as the range of / in a 
two-element sample of /, the first element chosen from the 
parent population of / in 2002-3, and the second element 
chosen from the same population in 2003-4. For a normal 
distribution, an estimate of the standard deviation can be 
obtained by multiplying the range by a coefficient k which 
is a function of the number of elements in the sam ple, n. 
Forn = 2, Table 12 of lCrow. Davis fc Maxfieldl Jl96Cll reads 
k — 0.886. Multiplying A/ by this value we get an estimate 
of the standard deviation of /. The latter number divided 
by the formal standard deviation of / from Table 2 yields 
the factor we are seeking. The mean value of the factor for 
the four frequencies turns out to be D = 2.0. If we applied 
the procedure to all nine frequencies, the result would be 
D — 1.8. If we used the formal standard deviations of the 
2002-3 solution, the results would be very nearly the same. 
Henceforth we shall adopt D — 2. A standard deviation 
equal to the formal standard deviation times this factor we 
shall refer to as "corrected," and from now on we shall drop 
the adjective "formal," so that by "standard deviation" we 
shall mean "formal standard deviation." 

Multiplying the standard deviations of the 2003-4 am- 
plitudes by two, we find that the ratios of the amplitudes 
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Table 3. A comparison of the frequencies and amplitudes of the nine independent 
terms common to Table 2 of Paper I and Table 2 of the present paper. The 
differences, A, are in the sense '2003-4 minus 2002-3.' 



ID 


A/ [d- 1 ] 


AA U [mmag] 


AA V 


[mmag] 


AA y [mmag] 


h 


-0.000014 


± 0.000017 


-1.2 ± 0.27 


-0.4 


± 


0.20 


-0.2 ± 


0.20 


fa 


-0.000033 


± 0.000028 


0.7 ± 0.27 


0.8 


± 


0.20 


0.3 ± 


0.20 


h 


-0.000081 


± 0.000030 


0.9 ± 0.27 


1.0 


± 


0.20 


1.3 ± 


0.20 


h 


0.000055 


± 0.000035 


-1.2 ± 0.27 


-0.6 


± 


0.20 


-0.6 ± 


0.20 


ft 


0.00079 


± 0.00032 


-0.6 ± 0.27 


-0.4 


± 


0.20 


-0.5 ± 


0.20 


h 


0.00060 


± 0.00048 


-1.6 ± 0.27 


-0.9 


± 


0.20 


-1.0 ± 


0.20 


h 


0.00020 


± 0.00035 


0.2 ± 0.27 


0.3 


± 


0.20 


0.4 ± 


0.20 


h 


0.00066 


± 0.00071 


0.1 ± 0.27 


0.1 


± 


0.20 


0.0 ± 


0.20 


/a 


0.00039 


± 0.00040 


-2.3 ± 0.27 


-1.2 


± 


0.20 


-1.5 ± 


0.20 



to the corrected standard deviations become approximately 
equal to the signal-to-noise ratios defined in the previous 
subsection. This is a pleasant surprise, lending support to 
our value of D. (Using the numbers from Table 2, the reader 
can verify the approximate equality of S/N and A v /(2a v ), 
where a v is the standard deviation of A v . Note that while 
for low frequencies A v /(2a v ) is about 30% larger than S/N, 
the approximate equality improves for high frequencies. This 
seems to be reasonable in view of the 1/f decrease of the 
noise level in the periodograms.) Additional support for 
D = 2 comes from the fact that a different line of rea- 
soning applied to m ultisite data similar to ours has led 
lHandler et alJ (|2000) to the same value. 

Let us now return to Table 3. As can be seen from the 
table, \AA\ are the largest for /a, fa, and fa (in this order). 
Multiplying the standard deviations of AA by D = 2 we 
find that (1) the amplitude of the low frequency term has 
decreased by 4.3cr c in u, 3.0<r c in v, and 3.8a c in y, where by 
cr c we denote the corrected standard deviation of AA, (2) the 
amplitude of the i = 6 term has decreased by 3.0, 2.2, and 
2.5<7 C in it, v, and y, respectively, and (3) the amplitude of 
the i = 3 term has increased by 1.7, 2.5, and 3.2(j c in u, 
v, and y, respectively. For the remaining six terms, \AA\ ^ 
2.2a c . We conclude that the decrease of the amplitude of 
the i = A term is real, that of the i = 6 term may be real, 
while the amplitude increase of the i = 3 term is probably 
spurious. In the remaining cases, there is little or no evidence 
for amplitude variation. 

The five independent frequencies which appear in the 
present Table 2 but not in Table 2 of Paper I are (in the 
order of decreasing v amplitude) /b, fa, fio, fa, fy The low 
frequency /b is very nearly equal to the frequency /21, one 
of three low frequencies derived in Paper II from the radial 
velocities of v Eri. However, in our 2003-4 power spectra 
we did not see the other two low frequencies of Paper II. 
Frequency fa was listed in Table 3 of Paper I as that of one 
of several "possible further signals." Frequencies fao and f y 
are new. 

Frequency fa was found in Paper I and tentatively as- 
cribed to £ Eri (see the Introduction). We shall demon- 
strate in Sect. 3.3 that this frequency and fa, the smallest- 
amplitude independent frequency in Table 2, are both due 
to £ Eri. 

In addition to nine independent terms, Table 2 of Paper 
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Figure 2. Low-frequency part of the power spectrum of the v 
differential magnitudes '/1 Eri — £ Eri' prewhitened with 0.615 
d — , 2/^ rb and 6/ orb . Arrows indicate peaks at frequencies equal 

^ forh> "^forh ' ^^orb ' ^^orb ^forh' 



I lists 14 combination terms. In the 2003-4 power spectra 
we did not see three of them, viz., fi + fs, /1 + /a + ft 
and /1 + 2/2. On the other hand, we detected differential 
combination-terms fi — fa, fi + f% — fa and f\ — fa which 
were not found in Paper I. 



3.3 The comparison stars 

The data used for analysis were the differential magnitudes 
'fi Eri — £ Eri.' Power spectra were computed independently 
from the u and v data. No power spectra were computed 
from the less numerous y data. 

In the first-run u and v power spectra, the highest peaks 
occurred at the same frequency of 0.615 d _1 . This frequency 
is close to f[ , the only one found for ji Eri in Paper I (see the 
Introduction). The second and third runs yielded frequencies 
of 0.272 and 0.815 d _1 , again the same for u and v. The first 
of these numbers is close to twice the orbital frequency of fi 
Eri, f' oA = 0.13589 d _1 , while the second, to six times this 
frequency. The low-frequency part of the power spectrum of 
the v differential magnitudes l /i Eri — £ Eri' prewhitened 
with 0.615 d _1 , 2f' orh and is shown in Fig. [5] In this 

figure, the arrows indicate peaks at frequencies equal to /o r b, 
3/orbi 4/o rb , 5/o rb and 8f^ lh . Peaks at these frequencies are 
also present in the u power-spectrum. 

The occurrence of so many harmonics of the orbital 
frequency implies that the data contain a strongly non- 
sinusoidal signal of this frequency. The first possibility that 
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Figure 3. The u (top), v (middle) and y (bottom) differential 
magnitudes 'fj, Eri — £ Eri' prewhitened with /j are shown as 
a function of orbital phase of fi Eri. Phase zero corresponds to 
H JD 2452800. 



comes to mind is an eclipse. Fig. |3] in which the differential 
magnitudes '/i Eri — £ Eri' prewhitened with 0.615 d _1 are 
plotted as a function of orbital phase, shows that /i Eri is 
indeed an eclipsing variable. 

Returning to frequency analysis of the comparison-stars 
data, we rejected observations falling in the orbital phase 
range from 0.4 to 0.54, i.e., the data affected by the eclipse, 
and recomputed the power spectra. In the first-run u and 
v power-spectra, the highest peaks occurred at the same 
frequency of 0.615 d _1 as before. The highest peaks in the 
second and third u run were at 0.701 and 0.813 d _1 , respec- 
tively, while the second and third v runs yielded 1.206 and 
0.701 d _1 . 

In the next step, we carried out nonlinear least-squares 
solutions separately for the u, v and y data. As starting 
values, we used all four frequencies found above, i.e., 0.615, 
0.701, 0.8132 and 1.206 d" 1 . The results are presented in 
the first four lines of Table 4. The fifth frequency, f 5 , will be 
explained shortly. The frequencies and their standard devi- 
ations, listed in column 2, were computed as straight means 
from the separate solutions for the three bands. The am- 
plitudes, A u , A v and A y , and their standard deviations, are 
given in columns 4, 5 and 6. The u-amplitude signal-to-noise 
ratio, computed in the same way as in Sect. 3.1, is listed in 
the last column. 

In the power spectra of the u and v differential mag- 
nitudes 'jj, Eri — £ Eri' prewhitened with the first four fre- 
quencies of Table 4, the highest peaks occur at the same 
frequency of 0.659 d _1 . Since a term of very nearly the same 
frequency is prominent in the frequency spectra of the 2002- 
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Figure 4. Power spectra of the u (lower panel) and v (upper) 
out-of-eclipse differential magnitudes 'fi Eri — £ Eri' prewhitened 
with the five frequencies of Table 4. Arrows indicate frequencies 
/ x = 10.874 and f y = 17.254 d _1 . 



3 comparison-stars data (see Sect. 4.2), we conclude that the 
signals at 0.659 d _1 are intrinsic. A five-frequency nonlinear 
least-squares solution yielded the value of the fifth frequency 
and the corresponding amplitudes given in the last line of 
Table 4. 

The power spectra of the u and v differential magni- 
tudes l n Eri — £ Eri' prewhitened with the five frequencies 
of Table 4 are shown in Fig. 2] In both spectra the highest 
peak occurs at 1.182 d _1 . Although these peaks may repre- 
sent another term in the variation of fj, Eri, we shall termi- 
nate the analysis at this stage for fear of over-interpreting 
the data. 

At high frequencies, peaks at / x = 10.874 and f y = 
17.254 d _1 can be clearly seen in the v spectrum (Fig. [I] 
upper panel). The -u-amplitude signal-to-noise ratio is equal 
to 5.2 and 4.4 for f x and f y , respectively. Although in the u 
spectrum in the lower panel the peaks at these frequencies 
are masked by noise, a closer examination shows that they 
are present. For / x , the it, v and y amplitudes computed 
from the differential magnitudes 'fj, Eri — £ Eri' amount to 
0.9 ± 0.22, 1.2 ± 0.16 and 1.0 ± 0.16 mmag (formal sigmas), 
respectively. To within one formal a, these numbers agree 
with the / x amplitudes derived from the l v Eri — £ Eri' 
differential magnitudes (see Table 2). For f y , the u, v and 
y amplitudes computed from the differential magnitudes 
Eri - £ Eri' are equal to 0.9±0.22, 0.9±0.16 and 0.7±0.16 
mmag, respectively. Again, there is a la agreement with the 
amplitudes obtained from the 'v Eri — £ Eri' data (see Table 
2). In addition, the phases agree to within la in all cases. 
We conclude that both frequencies are due to an intrinsic 
variation of £ Eri. 



4 LOW FREQUENCIES FROM THE 2002-3 
DATA 

4.1 The eclipse 

We have to admit that in our original analysis of the 2002-3 
data we missed the eclipse of Eri (see Paper I and the 
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Table 4. Frequencies, periods and amplitudes in the out-of-eclipse differential magnitudes l /i Eri — 
£ Eri' from the 2003-4 data. The last column contains the ^-amplitude signal-to-noise ratio. 



ID Frequency [d 



Period [d] 



[mmag] A v [mmag] 



[mmagj 



S/N 



0.61504 ± 0.00010 
0.70160 ± 0.00027 
f 3 0.81351 ± 0.00026 
/I 1.20739 ± 0.00027 



1.6259 ± 0.00026 
1.4253 ± 0.00055 
1.2292 ± 0.00039 
0.8282 ± 0.00019 



9.4 ± 

4.0 ± 
3.3 ± 

3.1 ± 



0.22 
0.22 
0.22 
0.23 



6.1 ± 
2.4 ± 

2.2 ± 

2.3 ± 



0.17 
0.16 
0.17 
0.17 



5.7 ± 0.16 

2.2 ± 0.15 
2.5 ± 0.17 

2.3 ± 0.16 



11.1 
3.9 
4.1 
6.0 



/' 0.65934 ± 0.00028 1.5167 ± 0.00064 3.4 ± 0.21 2.3 ± 0.16 2.4 ± 0.15 4.2 



Introduction). Fig. |K| shows phase diagrams in which the 
2002-3 differential magnitudes 'pi Eri — d; Eri' prewhitened 
with f[ are plotted as a function of orbital phase. The eclipse 
can be seen clearly. 

A comparison of the phase diagrams in Figs. [2] and [3] 
shows that while the middle of the eclipse in 2002-3 falls at a 
phase of about 0.30, in 2003-4 it does at about 0.47, indicat- 
ing a problem with Hill's (1969) value of the orbital period. 
Assuming that Hill's value yields a correct cycle count be- 
tween the first eclipse observed in 2002 and the last one in 
2003, we arrive at the following ephemeris: 



Min. light = HJD 2452574.04 (4) + E/0. 135490 (18), 



(4) 



where E is the number of cycles elapsed from the epoch 
given (which is that of the middle of the first eclipse we 
caught in 2002), and the numbers in parentheses are es- 
timated standard deviations with the leading zeroes sup- 
pressed. The question why our photometric period differs 
from Hill's spectrographic one will be addressed in a forth- 
coming paper. 



4.2 Analysis of the out-of-eclipse fi Eri data 

In 2002-3 the numbers of differential magnitudes 'fi Eri — £ 
Eri' in the three bands were nearly the same, amounting to 
2823, 2830 and 2919 in u, v and y, respectively. After we re- 
jected observations falling within the eclipse, these numbers 
were reduced to 2597, 2603 and 2688, still sufficient for anal- 
ysis. Using these reduced data, we computed power spectra 
separately for u, v and y. The first two runs yielded the same 
frequencies of 0.616 and 0.701 d -1 in all three bands. These 
frequencies are very nearly equal to f[ and of Table 4. In 
the third run, the highest peak in the u power-spectrum was 
at 0.657 d~ , while in the v and y power-spectra the highest 
peaks were at the same frequency of 1.207 d~\ The first 
of these numbers is close to f' 5 , while the second is nearly 
identical with fi (see Table 4). 

The fourth run was, however, a disappointment. In the 
u power-spectrum, the highest peak occurred at 1.000 d _1 , 
while in the v and y ones, at 0.032 d~ . Since neither of these 
frequencies is likely to be intrinsic, we did not attempt to 
compute fifth-run power spectra. 

The four frequencies found above, i.e., 0.616, 0.701, 
1.207 and 0.657 d _ were used as starting values in a four- 
frequency nonlinear least-squares solutions. Justification for 
including fi in the u solution comes from the fact that a peak 
at this frequency is prominent in the fourth-run u power- 
spectrum. Likewise, f'$ was included in the v and y solutions 
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Figure 5. The u (top), v (middle) and y (bottom) 2002-3 differ- 
ential magnitudes 'fi Eri — £ Eri' prewhitened with /J are shown 
as a function of orbital phase of /i Eri. As in Fig. |3] phase zero 
corresponds to HJD 2452800. 



because prominent peaks at this frequency can be seen in 
the fourth-run v and y power-spectra. 

The results of the four- frequency solutions are given in 
Table 5 above the horizontal line. This table has the same 
format as Table 4. However, the S/N (last column) was now 
computed from the y data. 

The four- frequency solutions did not include f$. Since 
peaks at this frequency were present in the fourth-run u and 
v power-spectra, we carried out a five-frequency nonlinear 
least-squares solutions for all five frequencies of Table 4. The 
resulting values of fg, the amplitudes and S/N are given in 
Table 5 below the horizontal line. 



4.3 v Eri 

As explained in Paper I, the 2002-3 differential magnitudes 
of v Eri were computed as l v Eri minus the mean of com- 
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Table 5. Frequencies, periods and amplitudes in the out-of-eclipse differential magnitudes '/i Eri — 
£ Eri' from the 2002-3 data. The last column contains the y-amplitude signal-to-noise ratio. 



ID Frequency [d 



Period [d] 



[mmagj 



mmag] A y [mmag] S/N 



f[ 0.61587 ± 0.00013 1.62372 ± 0.00034 9.9 ± 0.26 6.2 ± 

f 2 0.70143 ± 0.00021 1.42566 ± 0.00043 6.9 ± 0.25 4.4 ± 

f' A 1.20690 ± 0.00030 0.82857 ± 0.00021 3.2 ± 0.25 3.4 ± 

f h 0.65751 ± 0.00036 1.5209 ± 0.0008 4.2 ± 0.25 2.3 ± 

fL 0.8147 ± 0.0006 1.2275 ± 0.0009 3.1 ± 0.24 2.3 ± 0.19 0.8 ± 0.15 1.4 



0.19 4.9 ± 0.15 9.4 

0.19 3.3 ± 0.15 4.1 

0.19 2.7 ± 0.15 4.5 

0.19 2.2 ± 0.15 4.2 



parison stars,' but with the low- frequency variations of [i 
Eri filtered out. Thus, a peak at low frequency in the power 
spectrum of these differential magnitudes — if not caused by 
noise — would be due to an intrinsic variation of v or £ Eri. In 
the latter case, however, the peaks would be suppressed by 
a factor of four, while the corresponding amplitudes would 
be divided by two. 

In Paper I we found a prominent peak at 0.254 d _1 
in the amplitude spectrum of the v differential-magnitudes 
prewhitened with all 23 frequencies identified from the 2002- 
3 data. The power spectrum of the same data also shows 
a prominent peak at this frequency. However, there is lit- 
tle power at this frequency in the 2003-4 spectra shown in 
Fig. Q On the other hand, in the 2002-3 u, v and y power- 
spectra there are peaks at 0.615 d _ , a frequency very nearly 
equal to /b found in Sect. 3.1 from the 2003-4 data. The 
2002-3 it, v and y amplitudes at this frequency amount to 
2.5 ± 0.19, 1.2 ± 0.14 and 1.4 ± 0.12 mmag, respectively, in 
fair agreement with the 2003-4 amplitudes listed in Table 2. 
Since multiplying the 2002-3 amplitudes by two would make 
the agreement much worse, the possibility that /b is due to 
£ Eri can be rejected. 



5 ANALYSIS OF THE COMBINED DATA 
5.1 v Eri 

After slight mean-light-level adjustments, the 2002-3 and 
2003-4 differential magnitudes of v Eri were combined, sep- 
arately for u, v and y. The combined, 2002-4 data have the 
time base-line of 525.8 d. The analysis of the 2002-4 differ- 
ential magnitudes was carried out in the same way as that 
of the 2003-4 data (see Sect. 3). Sixteen independent and 20 
combination frequencies could be identified from the power 
spectra. In all cases but two the yearly aliases were signifi- 
cantly lower than the central peak, so that there was no ±1 
y _1 uncertainty. This was to be expected because the 2002-3 
and 2003-4 observing windows span as much as 0.43 y each 
(see Sect. 2). The two exceptions were fe and /12. They will 
be discussed later in this section. 

The 36 frequencies derived from the combined data are 
listed in the first column of Table 6. As in Tables 2, 4 and 5, 
the values of the independent frequencies and their standard 
deviations, given in column 2, were computed as straight 
means from the separate solutions for u, v and y. The com- 
bination frequencies, listed below the horizontal line, were 
computed from the independent frequencies according to ID 
in the first column; their standard deviations were computed 



from the standard deviations of the independent frequencies 
assuming rms propagation of errors. The amplitudes, A u , 
A v and A y , given together with their standard deviations in 
columns 3, 4 and 5, respectively, are from the independent 
solutions for u, v and y. The ^-amplitude S/N, computed in 
the same way as in Sect. 3.1, is given in the last column. It 
can be seen that all independent frequencies meet the signifi- 
cance condition of Paper I. Among combination frequencies, 
this condition is not satisfied in two cases, viz., fa + f<± and 

h - h- 

In addition to frequencies derived from the 2003-4 data 
(see Table 2), Table 6 contains two further high frequencies 
due to v Eri, viz., /n and /12. The latter is close to that of 
one of several "possible further signals" listed in Table 3 of 
Paper I and to frequency vr obtained in Paper III from ra- 
dial velocities of the Si III triplet around 457 nm. Frequency 
/n is new. In order to make sure that this frequency is not 
due to £ Eri, we examined the 2002-4 out-of-eclipse differ- 
ential magnitudes 'fi Eri — £ Eri' prewhitened with the six 
frequencies of Table 7 (see the next subsection). In the pe- 
riodograms, there were no peaks at /n; the highest peak in 
the vicinity, at 6.7168 d - , had the v amplitude equal to 
about 0.4 mmag and S/N < 2.5. Analogous tests with the 
2002-3 data also proved negative. 

We shall now discuss the two problematic frequencies, 
fs and /12, mentioned in the first paragraph of this section. 
In case of f%, the central peak at 6.2438 d _1 was only slightly 
higher than the +1 y _1 alias at 6.2465 d _1 . However, the 
2002-3 and 2003-4 values of fa are both much closer to the 
frequency of the central peak than to that of the alias peak. 
Moreover, in each band, the nonlinear least-squares fit con- 
verged to exactly the same solution regardless of whether the 
starting value of fe was the frequency of the central peak, 
the 2002-3 value, or the 2003-4 value. We conclude that f 6 
given in Table 6 is unlikely to be in error by 1 y _1 . 

The case of /12 is similar, but now the — 1 y _1 alias at 
6.2210 d _1 is the problem. In v and y it is almost as high as 
the central peak at 6.2236 d , while in u it is even slightly 
higher. Computing power spectra from the averaged u, v 
and y residuals, with proper weight given to each band, did 
not solve the problem. In Paper III, there is also the y _1 
uncertainty: the frequency is equal to 6.22304 d _1 for Si III 
455.3 nm, but for Si III 456.8 and 457.5 nm it is close to 
6.2210 d" 1 . The frequency given in Table 3 of Paper I is 
equal to the alias frequency. More data are needed to decide 
whether the value given in Table 6 is the correct one. 

Fig. |S| shows the power spectra of the 2002-4 data 
prewhitened with the 36 frequencies of Table 6. In the u and 
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Table 6. Frequencies and amplitudes in the differential magnitudes of v Eri from the combined 
2002-2004 data. Independent frequencies are listed in the upper part of the table. The combi- 
nation frequencies are listed below the horizontal line. In both cases the frequencies are ordered 
according to decreasing v amplitude, A v . Frequencies / x and f y are due to £ Eri. 



ID Frequency [d 1 ] A u [mmag] A v [mmag] A y [mmag] S/N 





h 


5.7632828 ± 


0.0000019 


72. 


.8 


± 


0. 


13 


10. 


.8 


± 


0. 


.10 


36. 


7 


± 


0. 


10 


214. 


7 




h 


5.6538767 ± 


0.0000030 


38. 


.5 


± 


0. 


13 


27. 


1 


± 


0. 


.10 


25. 


i 


± 


0. 


10 


142. 


6 




h 


5.6200186 ± 


0.0000031 


35. 


.0 


± 


0. 


13 


24. 


5 


± 





.10 


23. 


2 


± 


0. 


10 


129. 







fi 


5.6372470 ± 


0.0000038 


31. 


.8 


± 


0. 


13 


22. 


3 


± 





.10 


21. 





± 


0. 


10 


117. 


.4 




/-> 


7.898200 


± 


0.000032 


3 


.6 


± 


0. 


13 


2. 


6 


± 


0. 


.10 


2. 


.5 


± 


0. 


10 


11. 


5 




./a 


0.432786 


± 


0.000032 


4. 


.1 


± 


0. 


13 


2. 


.5 


± 


0. 


.10 


2. 


5 


± 


0. 


10 


8. 


3 




h 


6.262917 


± 


0.000044 


2. 


.8 


± 


0. 


14 


2. 





± 





.10 


1. 


8 


± 


0. 


10 


11. 







h 


6.243847 


± 


0.000042 


3 


.0 


± 


0. 


13 


1. 


9 


± 





.10 


2. 


1 


± 


0. 


10 


10. 


5 




fn 


0.61440 


± 


0.00005 


3 


.0 


± 


0. 


13 


1. 


.4 


± 


0. 


.10 


1. 


6 


± 


0. 


10 


5. 


5 




h 


7.91383 


± 


0.00008 


1 


.7 


± 


0. 


13 


1. 


1 


± 


0. 


.10 


1. 


2 


± 


0. 


10 


6. 


1 




/x 


10.87424 


± 


0.00012 





.8 


± 


0. 


13 


1. 





± 





.10 


0. 


7 


± 


0. 


10 


5. 


7 




/io 


7.92992 


± 


0.00010 


1 


.2 


± 


0. 


13 


0. 


9 


± 





.10 


0. 


9 


± 


0. 


10 


5. 







fs 


7.20090 


± 


0.00009 


1 


.1 


± 


0. 


13 


0. 


.9 


± 


0. 


.10 


0. 


.9 


± 


0. 


10 


5. 







in 


6.73223 


± 


0.00012 


1. 


.0 


± 


0. 


13 


0. 


.8 


± 


0. 


.10 


0. 


6 


± 


0. 


10 


1. 


5 




fi> 


6.22360 


± 


0.00012 





.9 


± 


0. 


13 


0. 


8 


± 





.10 


0. 


.8 


± 


0. 


10 


4. 


.4 




fy 


17.25241 


± 


0.00016 





.6 


± 


0. 


13 


0. 


6 


± 





.10 


0. 


5 


± 


0. 


10 


1. 


.4 


fl 


+ h 


11.4171595 


± 


0.0000036 


12. 


.4 


± 


0. 


13 


8. 


8 


± 





.10 


8. 


1 


± 


0. 


10 


50. 


9 


h 


+ /s 


11.3833014 


± 


0.0000036 


10 


.8 


± 


0. 


14 


7. 


6 


± 





.10 


7. 


3 


± 


0. 


10 


44. 





h 


+ /4 


11.4005298 


± 


0.0000042 


10 


.2 


± 


0. 


11 


7. 


2 


± 


0. 


.10 


6. 


.8 


± 


0. 


10 


11. 


7 




2/l 


11.5265656 


± 


0.0000027 


1 


.4 


± 


0. 


13 


3. 


1 


± 


0. 


.10 


2. 


9 


± 


0. 


10 


17. 


9 


+ h 


+ fs 


17.037178 


± 


0.000005 


3 


.6 


± 


0. 


13 


2. 


5 


± 





.10 


2. 


3 


± 


0. 


10 


18. 


1 


h 


+ fs 


11.2738953 


± 


0.0000043 


2. 


.6 


± 


0. 


13 


1. 


5 


± 





.10 


1. 


3 


± 


0. 


10 


8. 


7 


h 


-h 


0.1094061 


± 


0.0000036 


2. 


.6 


± 


0. 


13 


1. 


5 


± 





.10 


1. 


6 


± 


0. 


10 


4. 





2/i 


+ /2 


17.1804423 


± 


0.0000040 


1 


.9 


± 


0. 


13 


1. 


5 


± 


0. 


.10 


1. 


3 


± 


0. 


10 


10. 


9 


2/i 


+ /4 


17.163813 


± 


0.000005 


1. 


.7 


± 


0. 


13 


1. 


.4 


± 


0. 


.10 


1. 


2 


± 


0. 


10 


10. 


1 


2/i 


+ h 


17.1465842 


± 


0.0000041 


1 


.6 


± 


0. 


13 


1. 


3 


± 





.10 


1. 


1 


± 


0. 


10 


9. 


.4 


+ h 


+ h 


22.800461 


± 


0.000005 


1. 


.4 


± 


0. 


13 


1. 





± 


0. 


.09 


0. 


.8 


± 


0. 


10 


8. 


3 




2/2 


11.3077534 


± 


0.0000042 


0. 


.8 


± 


0. 


11 


0. 


.8 


± 


0. 


.10 


0. 


.5 


± 


0. 


10 


4. 


6 


+ h 


-h 


5.797141 


± 


0.000005 


1 


.0 


± 


0. 


13 


0. 


8 


± 





.10 


0. 


.8 


± 


0. 


10 


4. 


2 


/l 


+ h 


13.661483 


± 


0.000032 


1 


1 


± 


0. 


13 


0. 


7 


± 





.09 


0. 


.8 


± 


0. 


10 


4. 


5 


+ /3 


+ /4 


17.020548 


± 


0.000005 





.9 


± 


0. 


14 


0. 


.7 


± 


0. 


.10 


0. 


.7 


± 


0. 


10 


5. 


1 


h 


+ /4 


11.291124 


± 


0.000005 


0. 


.8 


± 


0. 


14 


0. 


.7 


± 


0. 


.10 


0. 


6 


± 


0. 


10 


1. 





+ h 


+ /4 


17.054406 


± 


0.000005 





.8 


± 


0. 


14 


0. 


6 


± 





.10 


0. 


8 


± 


0. 


10 


1. 


1 


h 


-h 


0.1260358 


± 


0.0000042 


1 


.7 


± 


0. 


13 


0. 


6 


± 





.10 


0. 


8 


± 


0. 


10 


1. 


6 


/i4 


-2/2 


17.071036 


± 


0.000005 





.6 


± 


0. 


14 


0. 


5 


± 





.10 


0. 


.4 


± 


0. 


10 


3. 


6 


h 


+ /4 


11.257266 


± 


0.000005 





.5 


± 


0. 


11 


0. 


.5 


± 


0. 


.10 


0. 


3 


± 


0. 


10 


2. 


9 



y spectra (bottom and top panel, respectively) the highest 
peaks occur at 0.3142 d~ , but in the v spectrum, at 0.2625 
d _1 . The corresponding amplitudes amount to 1.9 and 1.2 
mmag in u and y, and 1.1 mmag in v. In neither case does 
the signal-to-noise ratio exceed 3.4, so that these peaks are 
unlikely to be intrinsic. The reader may remember that in 
Sect. 3.1, a peak seen at 0.314 d _1 in the 2003-4 u power 
spectrum was also dismissed as spurious. 

The highest S/N peaks in Fig.HJhave frequencies equal 
to 5.0139 d- 1 in u (S/N = 4.5), 13.0152 d" 1 in v (S/N = 
4.2) and 22.9440 d" 1 in y (S/N = 4.1). The peak at 5.0139 
d _1 can be explained in terms of colour extinction in the u 
band (see Sect. 3.1). At 13.0152 d _1 , there are low peaks in 
the u and y spectra with S/N equal to 2.9 and 2.4, respec- 
tively. No peaks at 13.0152 d" 1 can be seen in the 'fj, Eri 
— £ Eri' power-spectra, so that this frequency is not due to 
£ Eri. Finally, the frequency of 22.9440 d _1 is very close to 
the combination frequency 3/i + fa . We conclude that while 



the latter frequency may be intrinsic, the former is probably 
spurious. 



5.2 Out-of-eclipse variation of Eri 

Combining the 2002-3 and 2003-4 out-of-eclipse differential 
magnitudes '/i Eri — £ Eri,' we obtained time series con- 
sisting of 5818, 5823 and 5189 data-points in u, v and y, 
respectively. 

The highest peaks in successive power-spectra of these 
data occurred at frequencies close to those found from the 
2003-4 and 2002-3 time-series separately (see Tables 5 and 
4) and at the frequency /g = 0.568 d _1 , which is new. In u, 
the frequencies appeared in the order f[, f^y fL f'i an( A /lai 
where the last frequency is the +1 y _1 alias of f' A — 1.2056 
d~ . The two highest peaks in the sixth-run power spectrum 
occurred at 2.009 and 0.997 d _1 , neither of which is likely 
to be intrinsic. The third peak, only slightly lower than the 
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Table 7. Frequencies, periods and amplitudes in the out-of-eclipse differential magnitudes l fi Eri 
from the combined, 2002-4 data. The last column contains the ^-amplitude signal-to-noise ratio. 



£ Eri' 



ID 



fx 

f'i 

Jl 
1% 



Frequency [d 



Period [d] 



A u [mmag] A v [mmag] A y [mmag] S/N 



0.615739 ± 0.000016 1.624065 ± 0.000042 9.9 ± 

0.700842 ± 0.000034 1.42686 ± 0.00007 4.8 ± 

1.205580 ± 0.000043 0.829476 ± 0.000030 3.0 ± 

1.208346 ± 0.000042 0.827578 ± 0.000029 2.9 ± 

0.658876 ± 0.000039 1.51774 ± 0.00009 4.0 ± 

0.56797 ± 0.00005 1.76066 ± 0.00016 2.7 ± 

0.81272 ± 0.00006 1.23044 ± 0.00009 2.8 ± 



0.16 


6.3 ± 


0.12 


5.6 ± 


0.11 


12.2 


0.16 


3.0 ± 


0.12 


2.5 ± 


0.10 


6.5 


0.16 


2.6 ± 


0.12 


2.4 ± 


0.10 


7.2 


0.16 


2.5 ± 


0.12 


2.4 ± 


0.10 


6.9 


0.16 


2.5 ± 


0.12 


2.6 ± 


0.11 


4.8 


0.16 


2.0 ± 


0.12 


2.1 ± 


0.11 


3.1 


0.16 


1.9 ± 


0.12 


1.4 ± 


0.10 


1.0 




Table 8. The mean separation, S, and asymmetry, A, of 
the close frequency triplets in the oscillation spectrum of v 
Eri. 



10 12 14 16 18 20 22 24 26 28 30 
Frequency [d 1 ] 

Figure 6. Power spectra of the combined, 2002-4 u (bottom), v 
(middle) and y (top) differential magnitudes of u Eri prewhitened 
with the 36 frequencies of Table 6. 



second one, was at /g. In v, the order was f[, f^, fi, fs, I'q 
and /3 , but in the third run the peak at was only slightly 
higher than the one at f' i& . In y, the order was the same as 
in v, except that in the last run the highest peak occurred at 
0.997 d~ , while the peak at f$ — although present — would 
be missed if it were not previously found in u and v. 

The 1 y _1 uncertainty which affects /i did not plague 
other frequencies; in all other cases the yearly aliases were 
significantly lower than the central peak. 

The results of the analysis are given in Table 7. The 
numbers for /i a are from nonlinear least-squares fits in which 
the alias frequency read off the power spectrum was used as 
the starting value. In these fits, the other frequencies and the 
corresponding amplitudes were only slightly different from 
those given in the table. 

The u to y amplitude ratio in Table 7 exceeds 1.2 for all 
frequencies. For four frequencies the ratio is greater than 1.5, 
while for and /g, it is equal to about 1.3. However, this 



Terms 


S [d- 1 ] 


A [d- 1 ] 


3,4,2 


0.0169290±0.0000022 


-0.000600±0.0000087 


12,6,7 


0.019658 ±0.000064 


-0.00118 ±0.00015 


5,9,10 


0.015860 ±0.00008 


0.00046 ±0.00019 



dichotomy may be illusory because the (formal) standard 
deviation of the latter number amounts to about 0.10. 

The reader may have noticed that f[ ~ /b. Since nei- 
ther frequency can be due to £ Eri because the amplitudes 
and phases do not match (see also Sect. 4.3), this curious 
near-equality must be an accidental coincidence. 



6 SUMMARY AND CLUES FOR 
ASTEROSEISMOLOGY 

6.1 Independent high frequencies of v Eri 

Fig. [7| shows schematically the 11 independent high- 
frequency terms of v Eri derived from the combined, 2002-4 
data. Comparing this figure with Fig. 4 of Paper I one can 
see that two of the three high-frequency "possible further 
signals" of Paper I are now upgraded to the status of cer- 
tainty. This has already been mentioned in Sect. 3.2 and 5.1. 
The terms in question are the i — 9 and 12 ones. Both are 
members of close frequency triplets. 

The third high-frequency "possible further signal" of 
Paper I, with frequency equal to 7.252 d , must remain 
in limbo. Although in the power spectra prewhitened with 
the 36 frequencies of Table 6 (see Fig. |SJ there is a series 
of low peaks in the vicinity of 7.25 d~ , the corresponding 
^-amplitudes are smaller than 0.6 mmag and the signal-to- 
noise ratios do not exceed 3.6. In u, the amplitudes are 
smaller than 0.8 mmag and S/N < 3.1. More data are 
needed to decide whether any of these peaks is intrinsic. 

For the close frequency triplets seen in Fig.0 the mean 
separation, S, and the asymmetry, A, axe listed in Table 8. 

In Sect. 5.1 we have warned that /12 may be in error 
by 1 y" 1 . If this were indeed the case, the values of S and 
A given in Table 8 for the 12,6,7 triplet should be replaced 
by 0.020964 and -0.00379 d"\ respectively. Because of this 
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5.5 6.0 6.5 7.0 7.5 8.0 
Frequency [d -1 ] 



much larger than the standard deviation of /a, but also 
larger than the frequency resolution of the 2002-4 data by a 
factor of about 1.5. In fact, 3/i — 3/3 coincides with the —1 
y _1 alias of /a- Since the alias has about the same height 
as the +1 y -1 one, the combination frequency's amplitude 
must be below the detection threshold. Consequently, there 
is no longer any doubt that /a is an independent low fre- 
quency in the variation of v Eri. In addition, we have found 
another low frequency, /b (see Tables 2 and 6). Since only 
high-order g modes have frequencies that low, the sugges- 
tion put forward in Paper I that v Eri were both a f3 Cep 
variable and an SPB star is amply confirmed. 



Figure 7. Schematic ti-amplitude spectrum of v Eri from the 
combined, 2002-4 data: the 11 independent high-frequency terms, 
numbered as in Table 6. 




350 400 450 500 550 600 
Wavelength [nm] 



Figure 8. Observed uvy amplitude ratios for the small-amplitude 
members of the 7.91 d — 1 triplet (circles with error bars) compared 
with theoretical amplitude ratios for I 2. The observed ratios 
are from the 2002-4 amplitudes (see Table 6), the theoretical ones, 
from Fig. 5 of Paper III. The open circles are shifted slightly in 
wavelength to avoid overlap. 

uncertainty, and the suspicion of a long-term variation of 
the amplitude of the i = 6 term (Sect. 3.2), the triplet is not 
particularly suitable for asteroseismology at this stage. For- 
tunately, the other two triplets are well-behaved. There are 
no y _1 uncertainties, even for the lowest-amplitude term 
of the 5,9,10 triplet, and no signs of long-term amplitude 
variation. In addition, the 1 = 1 spherical harmonic identi- 
fication for all members of the large-amplitude triplet and 
the i = 5 member of the 5,9,10 one are secure (see Paper 
III or the Introduction). As can be seen from Fig. |5] the 
i = 9 and 10 members of the triplet have uvy amplitude 
ratios consistent with I equal to 1 or 2. Unfortunately, the 
standard deviations of the amplitude ratios, especially those 
of the smallest-amplitude member, are too large to fix I un- 
ambiguously. 

6.2 Independent low frequencies of v Eri 

As we mentioned in the Introduction, the frequency resolu- 
tion of the 2002-3 data was insufficient to reject the possibil- 
ity that /a is equal to the combination frequency 3/i — 3/3. 
As can be seen from Table 6, the difference between /a and 
3/i — 3/3 amounts to 0.0030 d _ . This number is not only 



6.3 n Eri 

To the single frequency f{ = 0.616 d _1 , derived from the 
differential magnitudes 'a Eri — £ Eri' in Paper I, we add 
five further ones (see Table 7). The values of the frequen- 
cies and the decrease of the uvy amplitudes with increasing 
wavelength (for at least four frequencies) indicate that we 
are seeing high-radial-order g modes. Thus, as already sug- 
gested in Paper I, the star is an SPB variable. Note that 
rotational modulation, the second hypothesis put forward 
in Paper I for explaining the /( term is now untenable be- 
cause it does not account for multiperiodicity. 

As can be seen from Table 7, the periods P3, P 2 and 
P[ axe equally spaced, with the spacing equal to ~ 0.20 d, 
while Pi precedes P3 by twice this value. The equal spac- 
ing in period may be a manifestation of the well-known 
asymptotic property of high-order g modes of the same I. 
There are, however, the following two problems with this 
idea: (1) the term half-way between the P3 and P4 ones 
is missing, and (2) the period-spacing is rather large. Bet- 
ter data may solve the first problem if the missing term is 
simply too weak to be detected in our data. The second 
problem requires a comparison with the theory. Unfortu- 
nately, the only SPB-star model available in the literature 
JPziembowski. Moskalik fc Pamvatnvkhlll993T) has M = 4 
Mq, logi/L/Q = 2.51 and X c = 0.37, whereas /i Eri is more 
massive (by ~ 2 Mq), more luminous (by ~ 0.8 dex), and 
more evolved (see Paper I) . In the model, the largest period 
spacing (for £ = 1) is equal to ~ 0.07 d, almost a factor 3 
smaller than observed in \i Eri. Whether this disagreement 
can be alleviated with a model which better matches the 
star remains to be verified. If this turns out to be unsuccess- 
ful, one can still invoke the unlikely idea that an unknown 
amplitude limitation mechanism is suppressing the modes 
halfway between the observed ones. 

The possibility that instead of an equally-spaced period 
triplet, P3, P2, P{, we have a rotationally-split frequency 
triplet, f[, f 2 , /3, is much less likely. Indeed, for an I = 1 g- 
m ode with frequency equal t o 0.7 d" 1 in the SPB-star model 
of lDziembowski et alJ (Il993t) . the observed mean separation 
of the frequency triplet, equal to 0.09849 ±0.00003 d~\ leads 
to equatorial velocity of rotation, v e , of about 30 kms - (see 
their Fig. 13), whereas available estimates of v D sini of [i Eri 
range from 150 to 190 kms -1 (see Paper I). Increasing the 
model's radius in order to better match [i Eri may increase v e 
to about 60 kms -1 , still much less than the observed values. 
An additional problem is posed by the large asymmetry of 
the frequency triplet. The asymmetry is equal to 0.0268 ± 
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0.000 1 d , while the rot ational splitting seen in Fig. 13 of 
iDziembowski et alJ dl993T) is nearly symmetric. 

In addition to confirming the SPB classification of Pa- 
per I, we have found fi Eri to be an eclipsing variable. As can 
be seen from Figs. [3] and the eclipse is a transit, probably 
total, the secondary is fainter than the primary by several 
magnitudes, and the system is widely detached. As far as 
we are aware, the only oth er eclipsing variable with similar 
properties is 16 (EN) Lac jJerzvkiewicdlT979l) . except that 
in the latter case the eclipse is partial. (Interestingly, the 
discovery of an eclipse of this well-known j3 Cephei variable 
was a by-product a three-site campaign undertaken for ob- 
serving the star's pulsations.) Solving the /i Eri system will 
yield the primary's mean density and its surface gravity. 
This, however, is beyond the scope of the present paper. 

6.4 £ Eri 

The frequencies f K = 10.8742 and f y = 17.2524 d" 1 (see 
Tables 2 and 6) and the MK type of A2 V (see the Intro- 
duction) leave no doubt that the star is a S Scuti vari- 
able. The Stromgren indices, c\ = 1.076 and b — y = 
0.038 faauck fe Mermilliodl Il993) . are not reddened. This 
is not inconsistent with the star's Hipparcos parallax of 
15.6 6 ± 0.80 mas. Using the par allax and the V magnitude 
from lHauck fe Mermilliodl <199Sh one gets M v = 1.12±0.11, 
a value which, together with the b — y index, places the star 
about 0.02 mag to the blue of the observational blue edge 
of the S S cuti instability strip in the Mv vs. b — y diagram 
(see, e.g., iHandleJ 12002'). Apart from indicating the need 
for a slight revision of the blue edge, this position in the 
diagram suggests marginal pulsation driving as a possible 
explanation for the small uvy amplitu des. However, in view 
of the star's high v sin i of 165 kms" 1 iAbt fe Morreljl99Eft . 
anothe r expla nation may be provided by the hypothesis of 
iBreeerl il982Tl that fast rotation is a factor in limiting pul- 
sation amplitudes. 

Unfortunately, with only two small-amplitude modes 
the asteroseismic potential of £ Eri is insignificant. 
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